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Chapter  1 
INTRODUCTION 


"Saart  Boabs”  or,  aore  properly.  Precision  Guided  Munitions  (PGMs)  were  first 
Introduced  during  the  war  In  Southeast  Asia.  Both  the  Electro-Optical  Guided  Boab 
(EOGB )  and  the  Laser  Guided  Boab  { LGB )  were  auch  aore  successful  at  destroying  eneay 
targets  than  conventional  unguided  Iron  boabs.  Since  then,  nuaerous  other  PGMs  have 
been  designed,  developed,  and  deployed,  each  eaploylng  a  particular  segaent  of  the 
electroaagnetlc  spectrua  for  guidance.  Because  successful  use  of  these  weapons 
depends  on  a  wide  variety  of  ataospherlc  variables  (l.e.,  clouds,  haze,  lllualnatlon 
level,  etc.).  Air  Heather  Service  (AWS)  had  to  develop  new  weather-support 
techniques  for  these  weapons. 

This  report  deals  with  weather  support  to  a  particular  kind  of  PGM,  naaely, 
those  that  eaploy  a  TV  sensor  aounted  In  the  nose  of  the  weapon.  Specifically,  what 
follows  Is  an  evaluation  of  the  largely  statistical  techniques  developed  by  Mr  Ralph 
Huschke  of  RAND  Corporation  (Huschke,  1976).  While  the  Huschke  aethod  'permits  the 
direct  evaluation  of  weapons  systea  performance  In  any  weather  situation  or  cllaatlc 
reglae  that  can  be  extracted  as  a  subset  of  the  weather  data  base,”  It  Is  not 
designed  specifically  as  an  operational  weather  support  technique.  Rather,  It  Is 
Intended  as  a  aethod  which  will  test  the  use  of  a  statistical  *best  fit”  between  the 
sky-ground  ratio  and  available  weather  data  to  Infer  contrast  transalttance.  Thus, 
the  aethod  Is  only  as  good  as  the  relationships  that  Huschke  gives  In  his  report 
(Huschke,  1976,  page  11),  l.e..  Its  ability  to  duplicate  what  goes  on  In  the  real 
world.  Its  operational  weather  support  utility  as  a  forecast  tool  Is  yet  to  be 
proven.  The  reader  should  be  thoroughly  familiar  with  the  techniques  described  In 
Huschke's  report  prior  to  reading  what  follows. 


Chapter  2 

BACKGROUND  INFORMATION 


To  successfully  support  TV-guided  PGM  operations,  one  must  employ  far  more 
sophisticated  techniques  than  those  commonly  used  In  conventional  aviation 
meteorology.  Ambient  and  predicted  future  Illumination  levels,  target-to-background 
contrast  levels,  and  the  temporal  and  spatial  variability  of  atmospheric  haze  are  a 
few  of  the  factors  which  contribute  to  PGM  guidance  and  are  not  a  part  of  routine 
aviation  weather  support.  A  complete  discussion  of  the  complexities  of  such 
specialized  support  Is  available  In  the  AWS  Electro-Optical  (EO)  Handbook  (Cottrell, 
et.al.,  1979)  and  will  not  be  repeated  here. 

The  AWS  TV-Guided  Maverick  Missile  Weather  Support  Plan  (WSP)  tasks  the  USAF 
Environmental  Technical  Applications  Center  (USAFETAC)  to  "Provide  Phase  I  technique 
development  to  produce  the  appropriate  forecast  methods  and  determine  the  forecast 
skill  levels."  Since  the  adoption  of  this  WSP  and  with  the  support  of  the  EO 
Systems  Working  Group  (SWG)  USAFETAC  has  culled,  from  the  scientific  literature,  a 
variety  of  atmospheric  models  which  were  designed  to  simulate  the  loss  of  visible 
light  contrast  between  a  ground  target  and  an  airborne  TV  sensor.  These  models  fall 
into  two  general  categories:  radiative  transfer  models  which  employ  rigorous 
scattering  theory  to  compute  contrast  loss  for  specified  geometrical  and 
environmental  conditions  and  the  empirical,  so-called  “cookbook,"  models  or 
techniques  which  employ  grosser  assumptions  than  the  above  theoretical  models  and 
often  Involve  manual  methods.  Examples  of  the  first  category  are  the  Air  Force 
Geophysics  Laboratory's  { AFGL )  FLASH  Model  (Collins  and  Blattner,  1970),  the  Draper 
Lab's  DART  Model  (Whltnay  and  Malchow,  1977),  and  the  AWS  Haze  Model  (Breitling  and 
Plllpowskyj,  1970).  Examples  of  the  faster  empirical  techniques  are  those  due  to 
Huschke  (1976),  the  AWS  EO  Handbook  (Cottrell,  et.al.,  1979),  and  the  5th  Weather 
Wing  Manual  Weather  Support  Procedures  (5WW,  1977). 

Each  model  or  technique  has  obvious  advantages  and  drawbacks.  The  FLASH  Model 
offers  a  rigorous  and  exact  treatment  of  the  radiative  transfer  process,  but  could 
not  be  used  to  support  PGM  operations  In  anything  approaching  real  time  because  it 
requires  enormous  amounts  of  computer  time.  The  DART  Model  is  some  30  times  faster, 
but  Is  still  too  slow.  The  AWS  Haze  Model  is  computationally  fast,  but  some  of  its 
restrictive  assumptions  (flat  earth,  no  multiple  scattering,  etc.),  may  limit  its 
use  under  certain  geometrical  and  atmospheric  conditions  (low  sun  angles,  shallow 
dive  angles,  low  visibilities).  The  "cookbook"  methods  seem  to  offer  obvious 
advantages.  They  are  fast  and  do  not  require  a  large  computer,  high  forecaster 
skill  levels,  or  advanced  knowledge  of  the  radiative  transfer  process.  But  are  they 
"good  enough"  to  support  AF  operations?  To  answer  this  question,  one  must  determine 
how  good  Is  "good  enough."  This,  in  turn,  is  dependent  on  how  well  AWS  can  forecast 
the  variables  that  go  into  these  simpler  models.  Clearly,  a  technique  which  depends 
greatly  on  a  visibility  forecast,  for  example,  is  of  limited  use  if  the  visibility 
forecasting  skill  level  is  nil. 

The  techniques  developed  by  Huschke  (1976)  hold  the  possibility  of  answering 
USAFETAC's  search  for  an  "appropriate  forecast  method"  to  support  TV-guided  PGMs. 
H1$  method  has  the  advantage  of  employing  common  meteorological  observations  plus 
estimates  of  target  backgrounds  and  can  be  employed  by  either  an  Individual  remote 
forecaster  or  a  centralized  computer  facility.  This  report  attempts  to  verify  the 
validity  and  comment  on  the  possible  future  operational  utility  of  the  Huschke  model 
by  comparing  Its  output  with  that  of  three  theoretical  computer  models  under  similar 
meteorological  and  geometrical  conditions.  The  Intercomparison  of  the  three 
theoretical  models  Is,  of  course,  involved.  Since  most  of  the  5th  Weather  Wing 
techniques  duplicate  the  Huschke  method,  the  results  of  this  study  will  also  apply 
to  them.  Neither  the  "rapid"  nor  the  more  extensive  methods  described  in  the  AWS  EO 
Handbook  (Cottrell,  et.al.,  1979)  will  be  discussed  here. 

While  the  Huschke  model  and  the  theoretical  computer  models  share  some  basic 
elementary  assumptions  and  equations,  each  model  treats  the  vertical  structure  of 
the  atmosphere  and  other  input  variables  in  a  different  manner.  These  differences 
complicate  the  comparison  process.  The  comparing  "apples  and  oranges"  aspects  of 
this  study  have  been  reduced  as  much  as  possible,  but  some  differences  are  endemic. 
A  discussion  of  these  differences  follows. 
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HORIZONTAL  VISIBILITY  (km) 


The  principal  feature  of  the  Huschke  model  (and  one  easy  to  mi s 1 nterpret )  is  the 
determination  of  the  vertical  structure  of  the  atmosphere  from  the  decision  tree 
(Huschke,  197b,  page  21).  Basically,  Huschke  views  the  atmosphere  as  being  composed 
of  two  layers:  a  lower  layer  with  the  surface  visibility  extending  to  the  top  of 
the  mixing  layer,  and  an  upper  layer  extending  to  the  top  of  the  atmosphere  with  a 
visibility  of  80  km.  The  top  of  the  mixing  layer  is  either  1500  meters,  200  meters 
for  clear  sky,  low  wind,  and  low  sun  conditions,  or  at  the  base  height  of  any  clouds 
between  200  meters  and  2500  meters.  Since  radiative  transfer  models  such  as  FLASH, 
DART,  etc.,  cannot  treat  partly  cloudy  conditions,  this  study  must  be  limited  to 
clear  skies  only.  For  a  TV  sensor  above  the  mixing  layer,  the  contrast 
transmittance,  Tc,  must  be  computed  for  the  slant  path  of  interest  for  each  layer, 
and  their  product  is  then  the  Tp  for  the  entire  path  length.  This  Huschke  “stair¬ 
step"  treatment  of  the  vertical  structure  of  the  atmosphere  is  different  from  all 
previous  runs  of  the  FLASH,  DART,  or  AWS  models,  which  assume  an  exponentially 
decreasing  extinction  coefficient  from  a  surface  value  based  on  the  surface 
visibility.  Both  a  lack  of  funds  and  the  nonavailability  of  a  CDC  6600  computer 
precluded  USAFtTAC's  modification  of  the  FLASH  software  to  treat  a  stair-step 
extinction  coefficient  profile.  Output  from  the  DART  model  for  specified  Inputs 
were  available,  but  the  lack  of  a  running,  debugged  version  of  the  program  prevented 
USAFETAC  from  modifying  the  program  to  include  stair-step  profiles.  Since  the  AWS 
model  was  easily  modified  and  fast  running,  both  stair-step  and  exponential  profiles 
were  used  to  compute  slant  path  contrast  transmittance  values  for  a  TV  sensor  at 
12,000  feet  AGL.  Figure  1  graphically  depicts  the  difference  between  a  stair-step 
(200-meter  or  1500-meter  mixing  depth)  profile  versus  an  exponentially  decreasing 
one  for  a  surface  visibility  of  23  km.  Note  that  for  a  mixing  depth  of  1500  meters 
the  visibilities  for  the  stair-step  profile  are  smaller  (extinction  coefficients 
larger)  than  the  exponential  profile  at  all  altitudes.  Thus,  one  would  expect  (a 
priori)  computed  Tc  values  to  be  smaller  using  the  stair-step,  1500-meter  mixing 
depth  profile.  This  subject  will  be  discussed  further  in  the  following  chapters. 

Another  bothersome  factor  in  trying  to  compare  the  Huschke  model  with  the  three 
radiative  transfer  models  is  how  each  treats  the  visible  spectrum.  While  the 
Huschke  model  computes  contrast  transmittance  values  that  purportedly  represent  the 
photopic  region  of  the  visible  spectrum,  the  other  models  compute  Tc  for  individual 
wavelength  intervals  in  the  visible  spectrum  (0.40  to  0.75  microns).  One  could 
compute  T  c  values  for  for  several  wavelength  intervals  and  then  numerically 
integrate  them  to  give  a  Tc  value  for  the  visible.  However,  this  would  involve 
unacceptable  computer  run  times  for  both  the  FLASH  and  DART  models.  A  far  more 
acceptable  solution  to  this  problem  is  to  compute  Tc  values  at  0.55  microns  and  then 
assume  that  such  values  approximate  to  an  acceptable  degree  those  for  the  entire 
photopic  region.  This  assumption  was  tested  using  the  AWS  model  to  minimize 
computer  run  times.  The  Tc  values  for  eight  wavelength  intervals  (0.40,  0.45,  0.50, 
0.55,  0.60,  0.65,  0.70,  and  0.75  microns)  were  computed  for  a  variety  of  geometrical 
inputs  and  background  reflectances.  These  values  were  then  numerically  intergrated 
(via  a  combination  of  Simpson's  and  Newton’s  3/8th  rules)  and  compared  with  the 
computed  values  at  0.55  microns.  Figure  2  shows  that,  within  the  assumptions  upon 
which  the  AWS  model  is  based,  computed  Tc  values  for  0.55  microns  do  closely 
reproduce  those  for  the  entire  visible  wavelength  region. 

The  azimuthal  angle  dependence  of  contrast  transmittance  is  another  factor  which 
is  treated  differently  by  the  Huschke  model  than  the  other  types  of  models.  The 
Huschke  model  computes  a  Tc  which  Is  not  a  function  of  the  sun-target-sensor 
azimuthal  angle  while  FLASH,  DART,  and  AWS  models  compute  a  Tc  which  is  a  function 
of  azimuthal  angle.  For  certain  sun  and  dive  angles,  this  difference  can  be 
considerable.  Figure  3  depicts  the  difference  between  Tc  values  for  azimuthal 
angles  of  0°,  90°,  and  180^  as  computed  from  the  FLASH  model  to  demonstrate  how 
large  such  differences  can  be.  Note  the  large  variation  of  Tc  values  with  dive 
angle  for  a  solar  zenith  angle  of  60°  and  with  dive  angle  for  a  solar  zenith  of  85°. 
To  facilitate  the  comparison  and  evaluation  of  the  models,  computed  values  of  Tc  for 
0°,  90°,  and  180°  solar  azimuthal  angles  for  FLASH,  DART,  and  AWS  were  averaged 
[((Tc(0°)  ♦  Tc(90°)  +  Tr  ( 180° ) )  /  3  ]  and  this  mean  value  was  then  compared  with  the 
Huschke-derived  values  for  the  same  geometry.  While  a  weighted  mean  would  probably 
have  had  more  validity,  no  acceptable  method  of  choosing  such  weights  could  be 
found. 
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Figure  3.  Variation  of  Contrast  Transmittance  with  Azimuth  Angle  for  FLASH  Model 


In  addition,  Huschke  computes  visible  extinction  coefficients  using  the  equation 
(i  *  3. 352/ VI slbll 1 ty ,  whereas  the  other  models  use  the  more  common  relationship  8  « 
3.912/vi slbl 1 Ity.  This  results  In  Huschke-1 nferred  visible  extinction  coefficients 
being  about  15  percent  lower  than  the  other  models  (see  Huschke;  1976,  pages  7  and 
8). 

The  other  Important  input  variables,  dive  or  depression  angle,  solar  elevation 
or  zenith  angle,  surface  albedo,  and  visibility  were  the  same  for  all  models. 
However,  Huschke  (1976)  treats  specific  values  of  these  variables  as  belonging  to  a 
particular  category  as  Input  to  an  algorithm  to  compute  the  sky-ground  ratio,  while 
the  other  models  employ  the  values  of  these  variables  directly  as  Input. 

In  summary,  this  study  will  compare  the  contrast  transmittance  values  as 
computed  from  the  FLASH,  DART,  and  AWS  models  for  clear  sky  conditions  and 

a.  One  altitude,  12,000  feet  AGL 

b.  Two  mixing  depths,  200  meters  and  1500  meters 

c.  Two  visibilities,  23  km  and  5  km 

d.  Three  solar  zenith  angles,  20°,  60°,  and  85° 

e.  Three  albedos,  0.06,  0.18,  and  0.80 

f.  Seven  dive  angles,  85°,  70°,  50<>,  30°,  20°,  10°,  and  7°  (0°  Is  level  flight, 
90°  is  straight  down). 

The  vertical  attenuation  coefficient  profiles  at  0.55  microns  for  23  km  and  5  km 
were  taken  from  Elterman  (1970  and  1968). 


Chapter  3 

CONTRAST  TRANSMITTANCE  COMPARISON  FOR  23  KM 


The  figures  which  follow  In  the  remainder  of  this  report  employ  a  common 
notation.  FL  denotes  the  results  of  the  Air  Force  Geophysics  Laboratory's  (AFGL) 
FLASH  Monte  Carlo  model  with  an  exponential  extinction  coefficient  profile  (Elterman 
1968).  OT  stands  for  the  Draper  Lab's  DART  mode)  with  the  same  profile  as  Input. 
HI  Is  output  from  the  Huschke  model  with  a  mixing  depth  of  200  meters  and  a  layered 
or  stair-step  extinction  coefficient  profile.  H2  represents  output  from  Huschke's 
model  for  a  mixing  depth  of  1500  meters.  A1  refers  to  the  results  of  the  AWS  Haze 
Model  with  an  exponential  (Elterman,  1968)  profile  while  A2  Is  for  a  stair-step 
extinction  coefficient  profile  with  a  mixing  depth  of  1500  meters. 

Figures  4,  5,  and  6  depict  the  variation  of  computed  contrast  transmittance  for 
the  various  models  for  an  albedo  of  0.18  and  solar  zenith  angles  of  20°,  60°,  and 
85°.  An  element  common  to  all  these  figures  Is  that  the  Huschke  model  values 
bracket  nearly  all  the  values  for  the  other  models  with  HI  showing  the  most 
optimistic  Tc  values  and  H2  showing  the  most  pessimistic.  This  trend  parallels  the 
total  optical  depth  below  the  sensor  for  the  various  models  with  HI  having  the 
smallest  and  H2  the  largest  values  of  all  the  models.  A  descrepancy  in  the  Huschke 
model  Is  apparent  in  Figures  5  and  6.  At  nearly  vertical  dive  angles  (85°)  the 
Huschke  model  predicts  Tc  values  which  are  considerably  less  than  for  a  70°  dive 
angle.  This  aspect  of  the  Huschke  model  Is  not  shared  by  any  other  model.  It  is 
due  to  the  way  that  F  delta  is  found  from  Tables  3  and  4  of  Huschke's  (1976)  report 
which  in  turn  are  based  on  measurements  of  the  sky-ground  ratio  published  by  Duff 
(  1972).  The  footnote  below  Table  4  restricts  F  delta  to  a  value  of  1.0  for  dive 
(depression)  angles  between  50°  and  82°.  Howeve.r,  for  dive  angles  greater  than  82°, 
F  delta  Is  allowed  to  take  on  values  larger  than  1.0.  This  laraer  value  of  F  delta 
can  result  in  a  larger  value  of  the  sky-ground  ratio  estimate,  ‘bb,  which  then  results 
la  the  calculation  of  a  Tc  value  for  85°  that  Is  less  than  that  for  70°  for  certain 
solar  zenith  angles.  While  Duff's  measurements  seem  to  support  this  anomaly,  all 
other  evidence  points  to  the  situation  whereby  the  contrast  transmittance  Is  at  a 
maximum  for  nadir  dive  angles  and  decreases  in  value  as  one  looks  toward  smaller 
dive  angles  near  the  horizon.  This  Is  one  aspect  of  the  Huschke  model  that  will 
require  further  investigation. 

Figures  7a  and  7b  more  clearly  show  the  differences  between  various  pairs  of 
models.  Plotted  are  the  percentage  differences 


Tcl  -  Tc2  x  100 
Tc1 

between  contrast  transmittance  values  for  each  of  several  dive  angles  for  a  solar 
zenith  angle  of  20°  and  an  albedo  of  0.18.  The  Huschke  model  with  a  mixing  depth  of 

200  meters  (HI)  agrees  best  with  the  AWS  model  (Al)  and  quite  well  with  the  FLASH 

and  DART  models,  with  a  gradual  fall  off  in  agreement  as  one  approaches  near-horizon 
dive  angles.  The  HI  shows  Tc  values  that  are  consistently  higher  than  FLASH  as  one 
would  expect  since  Its  surface-to-sensor  optical  depth  is  smaller.  On  the  other 

hand,  there  is  no  apparent  reason  why  Al  with  the  same  optical  depth  and  aerosol 
model  as  FL  should  show  consistently  higher  Tc  values  than  FLASH.  The  Huschke  model 
with  a  1500-meter  mixing  depth  ( H 2 )  compares  less  favorably  with  FLASH,  especially 
at  dive  angles  less  than  30°.  Once  again,  its  greater  optical  depth  accounts  for 
some  of  the  variation.  Figure  7a  shows  that  the  percentage  difference  between  H2 
and  FLASH  and  DART  gradually  Increases  with  decreasing  dive  angle,  becoming  greater 
than  100  percent  near  the  horizon.  Figure  7b  again  shows  the  Huschke/AWS  comparison 
(Hl.Al)  along-side  the  comparisons  of  the  other  models  with  one  another.  Most 

notable  Is  the  fact  that,  as  the  optical  depth  Increases  for  smaller  dive  angles, 
the  AWS  model  compares  less  and  less  favorably  with  FLASH  and  DART  (although  not  as 
badly  as  Huschke).  Percentage  difference  plots  for  solar  zenith  angles  of  60°  and 
85°  show  similar  trends. 

Figures  8,  9,  and  10  show  the  variation  of  Tc  with  dive  angle  for  the  various 
models  for  an  albedo  of  0.06  and  solar  zenith  angles  of  20°,  60°,  and  85°.  Figures 
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Fiqure  4.  Contrast  Transmittance  at  20  Degrees  Solar  Zenith,  Albedo  of  0.18,  and 
23-km  Surface  Visibility. 
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Figure  0.  Contrast  Transmittance  at  60  Degrees  Solar  Zenith,  Albedo  of  0.18,  and 
23-km  Surface  Visibility. 
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Figure  6.  Contrast  Transmittance  at  85  Degrees  Solar  Zenith,  Albedo  of  0.18,  and 
23-km  Surface  Visibility. 


12 


*  DIFF  CONTRAST  TRANSMITTANCE 

100  -80  -60  -40  -20  0  20  40 


Figure  7b.  Model  Comparisons  at  20  Degrees  Solar  Zenith,  Albedo  of  0.18,  and  23-)tin 
Surface  Visibility. 


14 


Figure  8.  Contrast  Transmittance  at  20  Degrees  Solar  Zenith,  Albedo  of  0.06,  and 
23-km  Surface  Visibility. 
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Figure  9.  Contrast  Transmittance  at  60  Degrees  Solar  Zenith,  Albedo  of  0.06,  and 
23-km  Surface  Visibility. 
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Figure  10.  Contrast  Transmittance  at  85  Degrees  Solar  Zenith,  Albedo  of  0.06,  and 
23-km  Surface  Visibility. 
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9  and  10  again  show  the  discrepancy  of  the  Huschke  model  for  near  vertical  dive 

angles.  Figures  11,  12,  and  13  repeat  the  process  but  for  a  background  albedo  of  \ 

0.80  (approximating  snow  cover).  Note  that  for  low  sun  conditions  (85°,  Figure  13)  • 

the  Huschke  model  (HI)  now  underestimates  Tc  values  while  It  tended  to  overestimate  { 

such  values  relative  to  the  other  models  for  higher  sun  angles  (20°  and  60°) 

previously.  Figure  14a  compares  the  percentage  differences  between  the  two  Huschke 

outputs  (HI  and  H2)  with  the  ether  models  for  a  different  solar  angle  (60°).  Once 

again,  the  AWS  (Al)  model  compares  best  with  the  Huschke  200-meter  model  (HI)  and 

the  H2.Fl  and  Hl.DT  model  comparisons  show  differences  of  less  than  30  percent. 

Also,  the  Huschke  1500-meter  model  (H2)  compares  less  favorably  with  FL  and  UT  for 

dive  angles  less  than  30°  showing  differences  of  over  100  percent  for  dive  angles 

less  than  20°.  The  fact  that  all  differences  are  negative  for  all  model  comparisons 

for  an  85°  dive  angle  (Figure  14a)  can  again  be  attributed  to  the  sky-ground  ratio 

problem  at  near  nadir  dive  angles  mentioned  earlier.  Figure  14b  shows  that  the  AWS  j 

(Al)  model  compares  well  with  the  FL  and  DT  models,  predicting,  for  all  dive  angles, 

T r  values  that  are  more  optimistic  than  the  other  models.  Figures  15a  and  15b  show  j 

similar  trends,  except  for  a  solar  zenith  angle  of  85°.  Figure  15b  shows  that,  for 
situations  where  the  sun  is  near  the  horizon,  the  DART  model  does  not  agree  1 

particularly  well  with  FLASH  even  though  both  are  compatible  spherical  earth  models. 

While  the  AWS  (Al)  model  predicted  larger  values  than  FLASH  for  sun  angles  of  20° 

and  60°,  It  now  predicts  Tc  values  considerably  smaller  than  FLASH.  This  may  be  due 

In  part  to  the  fact  that  the  AWS  Haze  Model  Is  a  flat-earth  radiative  transfer 

model.  However,  It  does  a  better  job  of  duplicating  FLASH  Tc  values  at  small  dive 

angles  than  at  large  ones,  which  Is  difficult  to  explain  from  either  a  theoretical 

or  geometrical  point  of  view.  Additional  percentage  difference  plots  for  other  1 

angles  and  albedos  will  not  be  Included  because  they  would  be  repetitious  and  add 

little  to  this  report.  J 

Figure  16  depicts  a  method  of  comparing  the  various  models  for  all  dive  angles 
simultaneously  as  the  mean  of  the  squared  differences  between  Tc  values.  The  model 
comparison  with  the  smallest  root-mean-square  error  ( RMSE )  value  represents  the  best 
overall  Intercomparison.  Note  that  as  the  solar  zenith  angle  increases,  the  RMSE 
values  are  generally  smaller  because  of  the  effect  of  smaller  differences  between 
smaller  values  of  contrast  transmittance.  Lines  connecting  similar  comparisons  for 
different  sun  angles  demonstrate  the  sun  angle  effect.  For  example,  the  Al.FL  RMSE 
value  decreases  from  20°  to  60°  and  then  Increases  for  an  85°  solar  angle,  failing, 
however,  to  match  the  general  downward  trend  in  RMSE  values  mentioned  above.  The 
Huschke  200-meter/FLASH  ( HI  ,FL )  RMSE  values  continue  to  Increase  as  the  solar  zenith 
angle  Increases.  While  the  FLASH/DART  RMSE  values  are  smallest  overall  for  both  20° 
and  60°,  the  RMSE  values  Increase  threefold  at  85°. 

Since  all  dive  angles  were  weighted  equally  in  the  RMSE  values  for  Figure  16  and 
since  most  PGM  operations  take  place  at  shallow  dive  angles,  the  RMSE  values  were 
recomputed  with  the  Tc  values  for  85°  and  70°  dive  angles  omitted.  Figure  17  shows 
these  new  RMSE  values  as  compared  to  the  values  when  all  seven  dive  angles  were 
included.  Note  that  while  most  model  comparisons  show  larger  RMSE  values  when  the 
two  largest  dive  angles  were  omitted  (solar  angle  20°),  three  comparisons  (H2,A2; 

Hl.Al;  and  FL.OT)  show  smaller  values  for  a  solar  zenith  angle  of  60°.  However, 
both  the  Huschke  (Hi)  and  AWS  (Al)  models  evidence  large  Increases  in  RMSE  values 
when  compared  to  FLASH  at  60°  leading  one  to  expect  that  these  two  models  share  a 
common  deficiency.  This  trend  is  reversed  at  an  85°  sun  angle  as  both  Huschke  and 
AWS  evidence  reduced  RMSE  values  when  the  Tc  difference  values  for  85°  and  70°  dive 
angles  are  not  Included.  Such  evidence  might  suggest  that  both  these  models  may  be 
biased  toward  both  near-horizon  sun  and  dive  angles. 

In  summary,  both  the  Huschke  and  AWS  models  do  an  adequate  job  of  approximating 
FLASH  contrast  transmittance  values  for  a  surface  visibility  of  23  km.  The  Inherent 
limitations  Imposed  by  comparing  models  with  different  vertical  extinction  profiles 
(and  therefore  optical  depths)  probably  account  for  some  of  the  comparison 
discrepancies.  Figure  18  compares  the  total  optical  depth  values  (0.55  microns) 
below  the  aircraft  altitude  (12,000  feet)  for  an  exponentially  decreasing  extinction 
coefficient  profile  (FLASH,  DART,  AWS)  with  the  Huschke  model  for  mixing  depths  of 
200  meters  and  1500  meters.  The  total  optical  depth  for  the  Huschke  200-meter  model 
(0.203)  Is  based  on  a  visibility  of  23  km  for  the  lowest  200  meters  and  a  visibility 
of  80  km  from  there  to  3.66  km.  The  total  depth  for  the  Huschke  1500-meter  model 
(0.361)  Is  similarly  based  on  a  visibility  of  23  km  up  to  1500  meters  and  80  km 
above  that.  Note  that  as  the  mixing  depth  is  increased,  the  optical  depth  Increases 
until  It  equals  the  value  for  FLASH,  DART,  and  AWS  (0.230)  Just  above  400  meters. 
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Figure  12.  Contrast  Transmittance  at  60  Degrees  Solar  Zenith,  Albedo  of  0.80,  and 
23-km  Surface  Visibility. 
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Figure  13.  Contrast  Transmittance  at  85  Degrees  Solar  Zenith,  Albedo  of  0.80,  and 
23-km  Surface  Visibility. 
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Figure  14a.  Model  Comparison  at  60  Degrees  Solar  Zenith,  Albedo  of  0.18,  and  23-km 
Surface  Visibility. 
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Figure  14b,  Model  Comparison  at  60  Degrees  Solar  Zenith,  Albedo  of  0.18,  and  23-km 
Surface  Visibility. 
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Figure  lr>a.  Model  Comparisons  at  85  Degrees  Solar  Zenith,  Albedo  of  0.18,  and  23-km 
Surface  visibility. 
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Figure  15b.  Model  Comparisons  at  85  Degrees  Solar  Zenith,  Albedo  of  0.18,  and  23-km 
Surface  Visibility. 
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Figure  17.  RMSE  for  All  Solar  Zenith  Angles  Versus  RMSE  without  85-  and  70-Degree 
Solar  Zenith  Angles. 
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Figure  18.  Total  Optical  Depth  Below  12,000-Ft  Viewing  Altitude  for  23-km  Surface 
Visibility. 
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Thus,  one  would  expect  Tr  values  for  HI  to  be  lore  optimistic  than  the  other  three 
models  because  of  the  smaller  optical  depth  involved,  and  that  H2  would  produce  more 
pessimistic  Tc  values  because  of  the  larger  optical  depth  {0.361  versus  0.230).  The 
tendency  of  some  AWS  model  output  to  be  more  optimistic  than  FLASH  at  high  sun 
angles  and  pessimistic  at  near-horizon  sun  angles  1$  a  cause  for  some  concern. 
Also,  the  manner  In  which  the  Huschke  model  treats  the  sky-ground  ratio  for  dive 
angles  between  82°  and  90°  should  be  reexamined.  Considering  the  fact  that  the 
Huschke  and  AWS  models  are  hundreds  of  times  faster  than  FLASH  certainly  argues  In 
favor  of  further  testing  of  these  models  at  other  surface  visibilities.  Such  a 
comparison  for  a  surface  visibility  of  5  km  follows  In  the  next  chapter. 
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Chapter  4 


CONTRAST  TRANSMITTANCE  COMPARISON  FOR  A 
5-KM  SURFACE  VISIBILITY 


While  comparisons  of  radiative  transfer  models  at  larger  surface  visibilities 
(particularly  23  km)  are  fairly  common,  such  a  visibility  Is  not  representative  of 
urban  areas  of  the  eastern  US  or  central  Europe  Much  smaller  surface  visibilities 
usually  prevail,  particularly  In  the  winter.  A  5-km  visibility  was  chosen  as  the 
other  visibility  for  this  study  because  It  Is  much  more  representative  of  central 
Europe  and  probably  also  represents  the  lower  limit  of  observed  visibilities  that 
will  still  allow  successful  TV-guided  PGM  operations. 

The  radiative  transfer  model  designators  used  in  Chapter  3  (HI,  H2,  Al,  A2,  FL, 
DT)  apply  again  In  this  chapter  with  the  exception  of  the  DART  model.  No  output 
data  for  this  model  for  a  visibility  of  5  km  are  available  at  this  time,  or  will  be 
In  the  foreseeable  future.  Once  again,  the  exponential  extinction  coefficients  were 
taken  from  Elterman  (1970).  All  other  input  variables  were  the  same  as  before. 
Since  computed  contrast  transmittance  values  for  oblique  dive  angles  at  5-km 
visibility  tend  to  be  very  small,  any  Tc  values  less  than  0.001  were  considered  as 
Insignificant  and  were  omitted. 

Figures  19,  20,  and  21  show  the  variation  of  contrast  transmittance  with  dive 
angle  for  the  five  models  for  a  solar  zenith  angle  of  20°  and  background  albedos  of 
0.06,  0.18,  and  0.80,  respectively.  As  they  did  for  a  visibility  of  23  km,  the  two 
Huschke  models  bracket  the  Tc  values  for  FLASH.  The  AWS  exponential  model  (Al) 
evidenced  considerably  more  variability  with  albedo  than  before.  Once  again,  the 
total  optical  depth  for  FLASH  lies  between  the  values  for  HI  and  H2.  This  point 
will  be  discussed  more  at  the  end  of  the  chapter.  Al  predicts  slightly  smaller  Tc 
values  than  FLASH  for  an  albedo  of  0.06,  much  larger  values  for  an  albedo  of  0.18, 
and  then  slightly  larger  values  for  an  albedo  of  0.80.  Neither  the  200-meter  nor 
the  1500-meter  mixing  depth  Huschke  models  approach  the  FLASH  values  for  all  three 
albedos.  Figure  22  shows  some  of  these  differences  more  clearly.  While  the 
percentage  differences  between  Huschke  200  meter  (HI)  and  FLASH  are  fairly 
consistent  (albeit  large)  throughout  the  range  of  background  albedos,  the  AWS/FLASH 
(Al.FL)  differences  vary  considerably  with  albedo,  Indicating  a  possible  Instability 
In  the  AWS  model.  The  fact  that  Al  matches  FL  exactly  at  a  10°  dive  angle  for 
albedos  of  0.06  and  0.80  is  not  significant  because  of  the  small  Tc  value  (0.002) 
Involved. 

Figures  23,  24,  and  25  show  the  variation  of  contrast  transmittance  with  dive 
angle  and  albedo  for  a  solar  zenith  angle  of  60°.  Once  again  the  two  Huschke  (HI 
and  H2 )  and  AWS  (Al  and  A2)  models  bracket  the  FLASH  model  values  with  the  AWS 
exponential  model  (Al)  values  showing  the  same  variability  with  albedo  as  previously 
when  the  solar  zenith  angle  was  20°.  Only  at  a  background  albedo  of  0.80  does  the 
AWS  exponential  extinction  coefficient  model  approximate  the  FLASH  values.  Figure 
26  shows  that  none  of  the  models  does  a  good  job  of  matching  FLASH  at  all  dive 
angles  and  albedos,  with  the  Huschke  200-meter  model  (HI)  showing  a  consistent  bias 
toward  optimistic  Tc  values.  In  general,  the  percentage  differences  between  the 
Huschke  and  AWS  models,  vice  FLASH,  Increase  with  decreasing  dive  angle,  but  the 
smallness  of  the  Tc  values  for  near-horizon  dive  angles  make  really  meaningful 
comparisons  suspect. 

Figures  27,  28,  and  29  repeat  the  comparison  process  for  a  solar  zenith  angle  of 
85°.  The  Huschke  and  AWS  models  again  bracket  FLASH  for  all  three  albedos  with  the 
Huschke  (HI)  values  at  near-nadir  dive  angles  again  different  because  of  the  way  the 
model  handles  the  sky-ground  ratio  estimate  for  these  angles.  As  it  did  for  a  solar 
zenith  angle  of  60°,  the  AWS  exponential  model  (Al)  again  is  less  than  FLASH  for  an 
albedo  of  0.06,  greater  than  FLASH  for  an  albedo  of  0.18,  and  very  nearly  matches 
FLASH  for  an  albedo  of  0.80.  Except  for  nearly  nadir  viewing,  most  T?  values  for 
near-horizon  sun  angles  are  so  small  that  It  is  doubtful  that  even  high  contrast 
targets  could  be  acquired  even  at  low  aircraft  altitudes.  Figures  30  and  31  show 
that  neither  the  Huschke  nor  AWS  models  do  a  good  job  of  duplicating  FLASH  contrast 
transmittance  values,  although  the  AWS  models  show  slightly  more  consistency  and 
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Figure  19.  Contrast  Transmittance  at  20  Degrees  Solar  Zenith,  Albedo  of 
5-km  Surface  Visibility. 
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Figure  20.  Contrast  Transmittance  at  20  Degrees  Solar  Zenith,  Albedo  of  0.18,  and 
5-km  Surface  Visibility. 
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Model  Comparisons  at  20  Degrees  Solar  Zenith  and  5-km  Surface  Visibility. 
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Figure  23.  Contrast  Transmittance  at  60  Degrees  Solar  Zenith,  Albedo  of  0.06,  and 
5-km  Surface  Visibility. 
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Figaro  24.  Contrast  Transmittance  at  60  Degrees  Zoiar  Zenith,  Albedo  of  0.18,  and 
5-km  Surface  Visibility. 
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Figure  26.  Model  Comparisons  at  60  Degrees  Solar  Zenith  and  5-km  Surface  Visibility. 
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27.  Contrast  Transmittance  at  85  Degree 
rface  Visibility. 
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Figure  28.  Contrast  Transmittance  at  85  Degrees  Solar  Zenith,  Albedo  of  0.18,  and 
5-km  Surface  Visibility. 
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Figure  29.  Contrast  Transmittance  at  85  Decrees  Solar  Zenith,  Albedo  of  0.80 
5-ka  Surface  Visibility. 
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less  variability  than  Huschke.  Also,  for  these  sun  and  dive  angles,  the  AWS 
1500-meter  mixing  layer  stair-step  model  { A2 )  does  not  compare  favorably  with  the 
Huschke  1500-meter  model  (H2).  Both  do  badly  when  compared  to  FLASH. 

Figure  32,  which  shows  the  mean  of  the  squared  differences  (RMSE)  In  Tc  values 
for  all  dive  angles  indicates  that,  overall,  the  Huschke  model  compares  well  with 
the  AMS  model  at  all  sun  angles  although  neither  Is  particularly  good  when  compared 
with  FLASH. 

In  summary,  a  comparison  of  the  output  from  the  Huschke  and  AWS  models  with  that 
from  FLASH  with  similar  Input  reveals  that,  for  a  visibility  of  5  km,  neither  model 
Is  without  fault  if  FLASH  Is  assumed  to  be  ground  truth.  The  Huschke  model  Is 
either  overly  optimistic  or  overly  pessimistic,  depending  on  which  mixing  depth  is 
used,  and  the  AWS  model  output  is  highly  variable,  depending  on  which  background 
albedo  Is  used.  In  general,  both  the  Huschke  HI  and  AWS  A1  models  tend  to  predict 
much  larger  values  which  are  double  those  for  FL  In  some  cases.  Both  H2  and  A2 
routinely  predict  smaller  values  than  FL. 

Figure  33  compares  the  values  of  the  total  optical  depth  {0.55  microns)  below 
the  viewing  altitude  for  the  Huschke  (stair-step)  versus  the  FLASH,  DART,  and  AWS  A1 
models  (exponential).  Note  that  the  total  optical  depth  for  the  exponential  models 
(0.780)  lies  about  midway  between  that  for  the  Huschke  200-meter  model  (0.325)  and 
the  Huschke  1500-meter  model  (1.28).  As  the  mixing  depth  Increases  above  200 
meters,  the  total  optical  depth  Increases  until  It  equals  that  for  the  exponential 
models  near  a  mixing  depth  of  800  meters.  One  would  thus  expect  that  HI  would 
predict  much  smaller  Tc  values  (vice  FLASH,  DART,  AWS  Al)  and  that  H2  would  predict 
much  larger  values  than  these  three  models. 
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Fiqurc  32.  RMSE  Model  Comparisons  for  Albedo  of  0.18  and  5-km  Surface  Visibility. 
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Figure  33.  Total  Optical  Depth  Below  12,000-Ft  Viewing  Altitude  for  5-km  Surface 
Visibility. 


Chapter  5 

CONCLUSIONS  AND  RECOMMENDATIONS 


Conclusions 


As  Mentioned  previously  at  the  end  of  Chapters  3  and  4,  the  Huschke  and  AWS 
models  do  a  fairly  acceptable  job  of  duplicating  FLASH  contrast  transmittance  values 
when  the  surface  visibility  is  good  (23  km),  but  show  larger  deviations  from  FLASH 
as  the  visibility  decreases  to  5  km.  Some  of  these  deviations  can  be  attributed  to 
differences  In  the  total  optical  depth  below  the  sensor,  i.e.,  differences  between 
exponentially  decreasing  and  stair-step  extinction  coefficient  profiles.  A  less 
significant  portion  may  be  due  to  the  fact  that  the  Huschke  model  contrast 
transmi ttance  values  are  Independent  of  the  sun-target-sensor  azimuthal  angle.  If 
one  computes  Huschke  Tc  values  for  mixing  depths  at  which  the  exponential  and 
stair-step  profiles  optical  depths  are  equal  (425  meters  for  23  km  and  800  meters 
for  5  km)  the  Huschke  values  are  closer  to  the  FLASH  values  then  those  computed 
previously  for  mixing  depths  of  200  meters  or  1500  meters.  In  some  cases,  the  Tc 
values  now  fall  within  the  range  of  FLASH  Tc  for  0°,  90°,  and  180°  azimuthal  angles. 

The  Importance  of  being  able  to  predict  the  mixing  depth  and  surface  visibility 
with  a  fair  degree  of  accuracy  in  the  support  of  PGM  operations  cannot  be 
overemphasized.  Even  perfect  predictions  of  these  two  variables,  used  as  input  to 
the  Huschke  or  AUS  model,  would  still  not  result  In  perfect  predictions  of  slant 
path  contrast  transmittance  (assuming  FLASH  is  ground  truth).  When  one  further 
considers  that  predictions  of  mixing  depth  and  visibility  are  relatively  new  and  of 
uncertain  accuracy,  the  selection  of  a  model  to  support  TV-guided  PGM  operations 
becomes  more  complex.  Clearly,  a  supposedly  exact  radiative  transfer  model  (FLASH) 
with  Imperfect  Input  data  will  produce  contrast  transmittance  data  of  diminished 
usefulness  to  the  tactical  decision  maker.  On  the  other  hand,  an  approximate  model 
(Huschke,  AUS)  will  produce  data  that  may  be  less  exact,  but  at  a  large  reduction  In 
computer  processing  time.  The  price  one  pays  for  using  a  more  sophisticated  model 
to  support  PGMs  may  therefore  be  ill-advised,  since  the  three  most  important  Input 
variables  (mixing  depth,  visibility,  albedo)  are  based  on  still  unproven  forecasting 
techniques.  Therefore,  for  the  reasons  listed  above,  a  simple  model  such  as 
Huschke's  appears  to  meet  the  criteria  of  being  “good  enough"  to  support  TV-guided 
PGM  operations. 

Recommendations 


Based  on  the  results  of  this  study  and  the  arguments  presented  above,  USAFETAC 
makes  the  following  recommendations  relative  to  the  selection  of  a  model  to  support 
TV-guided  PGMs: 

a.  That  the  Huschke  model  be  modified  for  dive  angles  greater  than  82°  to 
produce  values  of  contrast  transmittance  that  more  nearly  approximate  the  FLASH 
predictions. 

b.  That  this  modified  Huschke  model  be  adopted  by  AUS  as  the  standard  model  for 
both  field  and  centralized  facility  support  to  TV-guided  PGM  operations. 

c.  That  AUS  units  employing  the  Huschke  model  undertake  an  ongoing  evaluation 
to  Improve  the  overall  capabilities  of  the  model. 
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GLOSSARY 


AF 

Air  Force 

AFGL 

Air  Force  Geophysics  Laboratory 

AGL 

above  ground  level 

AWS 

Air  W eather  Service 

Al 

AU'S  ha/e  model  with  exponential  profile 

A2 

stair-step  extinction  coefficient  profile  with 
1500  meters 

mixing  depth  of 

DART 

Oraper  Lab's  radiative  transfer  model 

OT 

Draper  Lab's  DART  model 

EO 

Electro-Optical 

EOGB 

Electro-Optical  Guided  Bomb 

fL 

resuits  of  FLASH  Monte  Carlo  model 

flash 

AFGL  radiative  transfer  model 

hi 

output  from  Huschke  model  and  mixing  depth  of 

200  meters 

H2 

output  from  Huschke  model  and  mixing  depth  of 

1500  meters 

LGB 

Laser  Guided  Bomb 

PGM 

Precision  Guided  Munitions 

SWG 

Systems  Working  Group 

Tc 

contrast  transmittance 

USAFETAC 

US  Air  Force  Environmental  Technical  Applications  Center 

WSP 

Weather  Support  Plan 

5WW 

5th  Weather  Wing 
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